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Abstract 



The construction of positronium decay amplitudes is handled through the 
use of dispersion relations. In this way, emphasis is put on basic QED princi- 
ples: gauge invariance and soft-photon limits (analyticity). 

A firm grounding is given to the factorization approaches, and some am- 
biguities in the spin and energy structures of the positronium wavefunction 
are removed. Non-factorizable amplitudes are naturally introduced. Their 
dynamics is described, especially regarding the enforcement of gauge invari- 
ance and analyticity through delicate interferences. The important question 
of the completeness of the present theoretical predictions for the decay rates 
is then addressed. Indeed, some of those non-factorizable contributions are 
unaccounted for by NRQED analyses. However, it is shown that such new 
contributions are highly suppressed, being of C {o-^) ■ 

Finally, a particular effective form factor formalism is constructed for para- 
positronium, allowing a thorough analysis of binding energy effects and ana- 
lyticity implementation. 



PACS Nos : 36.10.Dr, 12.20.Ds, 11. 10. St, 11.55.Fv 
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1 Introduction 



Positronium is a bound state made of an electron with a positron. In this paper, 
we will be interested in the ground states, parapositronium {p-Ps, singlet, J — 0) 
and orthopositronium (o-Ps, triplet, J = 1). Positronium is unstable, it decays 
when its components annihilate into photons. The lifetimes are quite different for 
the two species, because parity and charge conjugation conservations imply that 
p-Ps decays into an even number of photons, while o-Ps can only decay into an 
odd number. 

Positronium was discovered in 1951 [1]. Since then, many measurements of the 
lifetime of both species have been achieved. The most precise to date are 

T'^'P {p-Ps 27) = 7.9909 (17) x 10^ sec"^ [2] (1) 

( 7.0398(29) ^sec-i Tokyo (SiOa Powder) [3] 
{o-Ps ^ Sj) ^ I 7.0514(14) ^sec-i Ann Arbor (Gas) [| (2) 
[ 7.0482 (16) /.tsec^i Ann Arbor (Vacuum) g 

Obviously, the situation is not clear for the o-Ps lifetime, since results are mutually 
exclusive. This is known as the orthopositronium lifetime puzzle. 

The corresponding lowest order theoretical predictions were first obtained as 

^p-Ps = = 8.033 X lO^sec-^ [6] (3a) 

2 (tt^ - 9) , ^ 

To-Ps = a'^m^- ^ = 7.211 /isec'i § (3b) 

97r 

with a the fine structure constant and m the electron mass. A great deal of work 
has been done on the calculation of radiative corrections, which include perturbative 
QED corrections to a non-relativistic treatment of the bound state wavefunction. 
The present highly accurate theoretical predictions for the respective lifetimes are 
(see for example [§| to [^): 

Tp_ps = (1 + 62^ + 5i^) = 7.989620 (13) x 10^ sec'^ (4a) 

2 (tt^ - 9) 

To-Ps = a^m^- ^ (1 + 53^ + (^5^) = 7.039965(10) At sec"^ (4b) 

The corrections to the dominant 2j and 37 modes are given by 

62-y = -Ap- + 2a^\n- + Bp (-] - — In^ - + C„— In - 
TT a \TT/ In a tt a 

w -0.0058828 + 0.0005242 + 0.0000277 - 0.0000045 - 0.0000048 

(537 = -\n-+Bo(-] -— In - + Co— In- 

TT 6 a \iT/ Ztt a TT a 

w -0.0238939 - 0.0000873 + 0.0002402 - 0.0000045 + 0.0000034 
with coefficients 

Ap = 5- 7rV4 Ao = 10.286606 (10) 

Sp = 5.14(30) Bo = 44.52 (26) 

Cp = -7.919(1) Co = 5.517(1) 

The non-logarithmic O (a^) corrections Bp p and logarithmic O {a^ In a) corrections 
Cp^o have been obtained only recentl y [[l2[ . The small contributions to the lifetimes 
from multi-photon decay modes are ||l6|l 

Si-y = 0.274 (1 

^57 = 0.19(1) 



TT 

a 
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The four-photon mode increases the p-Ps width by about 1 x lO^sec"^, while the 
five-photon mode contributes for 7 x 10~^ /xsec~^ to the o-Ps width. 

As can be observed, agreement between theory and experiment is good for para- 
positronium, but the Bp, Cp radiative corrections are not yet accessible experimen- 
tally. For orthopositronium, these corrections render the theoretical prediction still 
closer to the experimental measurement of Ref. |^ (anyway, beyond discriminat- 
ing between experimental results, Bq and Co radiative corrections are not tested 
either). 

Positronium is a test ground for bound state treatments in Quantum Field The- 
ory. The first try dates back to the 40 's, with decay rates expressed through a 
factorized formula 

F {p-Ps ^ 27) = • {^vreia (e+e- ^ 27))^^^^^^ (5) 

with 0o the Schrodinger positronium fundamental state wavefunction at the origin, 
a (e+e~ 27) the total cross section for e+e~ 27 and Vrei the relative velocity 
of and in their center-of-mass frame. Since then, more sophisticated decay 
amplitudes have been constructed, and systematic procedures for calculating cor- 
rections have been developed (like non-relativistic QED (NRQED), see for example 
[17]), but factorization remains central. To avoid misunderstanding, let us stress 
that by factorization is usually meant a positronium decay amplitude built as a 
convolution integral, 

M (p-Ps ^ 27) - J d^k ij{k)-M (e+ (-k) e- (k) ^ 27) (6) 

with ip (k) the momentum Schrodinger wavefunction and M (e+e^ 27) the scat- 
tering amplitude for on-shell electron-positron into photons. 

Several questions concerning this formalism need to be addressed: 

- The basic factorization of the bound state dynamics from the anni- 
hilation process remains as the basic postulate. For low order correc- 
tions, this approximation is unquestionable, but for O (a^) corrections, 
factorization has to be tested. Indeed, non-perturbative phenomena 
responsible for the off-shellness of the electron and positron inside the 
positronium are of O (a^) (since positronium mass minus twice the elec- 
tron mass is of that order). In other words, to get a sensible theoretical 
prediction at O (a^), one must carefully analyze how binding energy 
effects enter the general factorization approach. 

- Violation of energy conservation ? Positronium being a bound state, 
M < 2m {M the positronium mass and m the electron mass). The 
total energy of the on-shell electron-positron pair entering the scattering 
process is 2vl? + n?, clearly greater than M, and getting worse as k 
increases in the integration. 

- The enforcement of gauge invariance may be problematic. This is 
more technical. Let us say that, in general, one has to project the 
electron-positron pair into a spin state compatible with the total spin 
of the bound state. This bound state spin state is usually treated non- 
relativistically, introducing some potential threads to the gauge invari- 
ance of the decay amplitude. 

- Some further problems concern the Bethe-Salpeter wavefunction. There 
are a lot of technicalities and subtleties associated with the bound state 
wavefunction. This wavefunction can in principle be obtained covari- 
antly from the Bethe-Salpeter (BS) equation. Usually, one has to in- 
troduce some approximations to solve this equation, ending with a non- 
covariant four-dimensional wavefunction. This 4Z) wavefunction is then 
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reduced to 3-D by making additional approximations for its temporal 
part (energy-dependent part). This last reduction may be questionable, 
(and, as said above, non-covariant wavefunctions can be difficult to ac- 
commodate with gauge invariance). 

- Last but not least, analytical behavior in the soft photon limit may be 
ill-defined. By this we mean that since factorization treats intermediate 
charged states as on-shell, infrared singularities are introduced. This 
problem is well-known, but some confusion seems to creep in the litera- 
ture. We stress that: For a model to be correct, it is not sufficient 
to show that IR divergences are successfully disposed of. From 
general quantum field theory principles, a decay amplitude involving 
only neutral external bosons has a correct analytical behaviour if and 
only if this amplitude vanishes in the soft photon limit (not simply be 
finite). This is a very profound principle, at the root of electrodynamics 




The purpose of the present paper is to answer those questions using dispersion 
relations (DR). As will become clear in the course of the presentation, DR appear 
as the natural framework to study positronium decay amplitudes. Some questions 
on the applicability of DR analyses to bound state decays could be raised, but, in 
our view, the fact that all the questions above find consistent answers gives great 
confidence in this approach. 

The most prominent achievements of DR analyses are 

1. a well-defined systematic procedure for factorizing the bound state 
dynamics from the annihilation process, 

2. an identification of the dynamics ensuring both gauge invariance and 
the correct soft photon limits. 

The basic idea of our approach is to consider the positronium decay amplitudes 
through a fully relativistic loop model, where only off-shell constituents appear. The 
nice surprise is that standard three-dimensional convolution-type amplitudes like (^ 
are recovered through a DR. Contrary to the usual approaches, no approximations 
are needed to reach the three-dimensional form. Most importantly, we will show that 
many variants of the standard formula involve some unnecessary and dangerous (i.e. 
gauge-dependent) approximations in their treatment of spins, and we will remove 
them. Since our derivation relies on well-established techniques of quantum field 
theory, we conclude that some O (a^) corrections have been forgotten in those 
works. 

It is when analytical properties are analyzed that our procedure gains its full 
respectability. We will see that DR introduce some additional non-factorizable (i.e. 
non-convolution-type) contributions to the decay amplitudes. Further, some bound 
state structure-dependent non-factorizable contributions must be introduced. All 
those additional contributions interfere with the standard amplitude contributions 
to enforce both gauge invariance and the correct analytical properties. Interestingly, 
because the present approach relies on a dynamical characterization of the various 
contributions, it permits an evaluation of their magnitudes. This is an important 
fact because if those additional contributions prove to have been missed in current 
calculations, the positronium lifetime O (a^) corrections would be wrong. We will 
not be able to provide a definite answer to this question, but nevertheless point 
towards an optimistic view. Provided no pathological enhancement occurs, the 
NRQED result indeed should contain all O (a^) corrections. 

The paper is articulated in three parts. The first part contains the presentation 
of the general results, and discussions about their implications. After a short review 
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of the standard approaches to the construction of positronium decay amphtudes, we 
go on to prove that our model reproduces those standard decay amphtudes. This 
result is completely general. We then apply the insight gained to the analysis of soft- 
photon behavior, and find the mechanism at play to ensure the correct analytical 
limits. 

The remaining of the paper contains applications and illustrations. We analyze 
in some details parapositronium decay into two photons in the second part of the 
paper. This decay is a bit special because the kinematics forbids soft photon energies 
to vanish. Also, the appearance of a pseudoscalar coupling (parapositronium is 
a pseudoscalar) has interesting consequences, allowing form factor-like dispersion 
relations to be built. This technique is then used in various calculations. This 
section closes by the description of an improved basis for the perturbation series 



(4a), where binding energy effects are singled out. 

The final part contains discussions about soft-photon limits. We take as an 
example the decay of paradimuonium (the /i^/i^ bound state, p-Dm) into e^e~7. 
This is the simplest positronium-like decay process where a photon can have a 
vanishing energy. Again, the pseudoscalar property of paradimuonium simplifies 
the discussion, allowing a thorough analysis of the non-factorizable contributions. 
We also briefly discuss the decay Ks — » e+e^7, modeled by an intermediate charged 
pion loop. The reason is the striking similarity between this "elementary" particle 
decay and the bound state analogue p-Dm — > e~^e~"f. For example, exactly the same 
mechanism is responsible for the implementation of the correct soft-photon limits 
(SPL). It is very interesting to see Quantum Field Theory at play in transforming 
soft-photon singularities from intermediate charged states into a vanishing SPL, no 
matter the details of the specific models. Also, we hope that this incursion into the 
realm of kaon physics will convince the most sceptical readers that our approach is 
viable and correct, as it is based on general quantum field theory principles. 

NB : This paper has emerged from three previous works, Ref. p9[ . 

2 Positronium Decay Amplitudes 
and Dispersion Relations 

Standard amplitudes for the decay of parapositronium (J = 0) or orthopositronium 
(J = 1) are defined as (see for example 0, jl^, [^ ) 



xM (e- (k,C) ,e+ {-k,a -> (7) 

where M is the bound state mass and ■0 (k) is the Schrodinger bound state wave- 
function. This wavefunction is taken in a convolution integral with the scattering 
amplitude for on-shell e+e^ A, constrained to the required spin state. By work- 
ing out this constraint, we reach the commonly quoted positronium decay amplitude 
(hereafter denoted 'standard approach' amplitude) 

M C'^'S -^A)=VWlj ^ V (k) ^Tr I i±j!pr (e- (k) , e+ (-k) ^ A) } 

(8) 

where 

M (e- (k, , e+ (-k, C) -^A)=v (-k, C) T {er (k) , e+ (-k) A) u (k, C) 
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P = 7^ for parapositronium, / for orthopositronium with polarization vector e^. 
Finally, the width is calculated as 

with 5 a symmetry factor for the final state A. 

This standard decay formula has many peculiar features, and needs a rigorous 
groimding. As quoted in the introduction, factorization and SPL properties arc 
the most pressing problems. In section 2.1, we intend to concentrate on the first 
point, leaving analytical questions to the following one. To do so we will particu- 
larize the discussion to parapositronium decay to 27, which cannot suffer from any 
pathological SPL due to its kinematics. 

Let us now introduce our model. Basically, we assume a loop structure for 
the decay amplitudes. In other words, positronium decays into a virtual electron- 
positron pair which subsequently annihilates into real or virtual photons (an odd 
number for ortho-states, an even number for para-states). The coupling of the 
positronium to its constituents is described by a form factor, denoted by Fb. It 
is not assumed to be a constant, since a constant form factor would amount to 
consider positronium as a point-like bound state. 

For parapositronium decay into two photons, our model is represented in figure 
1. The corresponding amplitude is 

>,(p-P5-77)= / -^FBTri^f,—^ ^""j-T^ j ^^^^2. (10) 

J (27r) I g-^r-m ^+2^-"^J 

with P the positronium four- momentum and Fb = Fb (^q^,P-q) the form fac- 
tor. The tensor F^'^ is the scattering amplitude for off-shell e^e^, with incoming 
momenta — q and |P -|- q, into two photons : 

= F"^ (e+ (iP - g) e- (iP + 5)^7(^1) 7 (M) 

= — r-^/ ieY + ieY .^ 1 ^ iej'' (11) 

RemEtrks : 

1) The model is extended to orthopositronium decays through the re- 
placement of 75 by ^. 

2 ) The electron and positron in the loop are never on-shell, because 
M < 2m. 

3 ) Fb contains all the information about the bound state. Let us 
postulate a form for this coupling as (in the positronium center-of-mass 
frame) 

Fb =CcPaT{qo,q^) {q^ +-f^) (12) 

with C a constant, and 7^ = — related to the binding energy 

Eb = 2m — M. In QED, Eb and 7^ are related to the fine structure 
constant as 7^ rs Tr?o? jA: and Eb = —ma^ /A. 

2.1 The loop model reproduces standeird decay amplitudes 

We can state our result in three steps (details are found in the appendix). 
First, we define 

Im T (p2) =lmM {p-Ps (P^) ^ 77) (13) 
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with the absorptive part given by the two vertical cuts, figure 2. 

Then, using an unsubstracted dispersion relation (see |2^) with s = P^: 

T(M').RoT{M^)^-£_^I„,TW (14) 

(T (M^) = ReT (M^) because M < 2m), and changing of variables, one recovers 
a factorized form (i.e. a convolution type amplitude): 

M{p-Ps^^j)^^ J -^^^:^|l^^rr{75(m-r)r(m+m (15) 

with r given by T'^'' (e" (fc) , e+ (fc') -> 77) e^^ej, and fc = (i;^, k) , fc' = (i;^, -k), 
i.e. the same F as in (8). 

Finally, if one further identifies C = ^/M/m and 0o^(O,k^) — ^(k), by ne- 
glecting the k dependence in the projectors (m— fi') and (m+ /(f), (15) reduces to 

®- 

As a corollary, the energy-dependence of the form factor is demonstrated to 
be irrelevant. Indeed, the dispersion relation fixes qo = 0. From now on, we will 
simplify the notation and identify T (O, k^) = !F (k-^) . 

2.1.1 Discussion 

Let us now comment on this result. Interpreting standard approach amplitudes as 
dispersion relation integrals, we can give answers to most of the questions raised in 
the introduction. 

First, factorization appears simply as a manifestation of the optical theorem: 
the appearance of on-shell intermediate states is expected in the imaginary part. 
The dispersive integral "shifts" them off-shell: the amplitude at the physical point 
s = JVP is purely real, i.e. only off-shell e+e" circulate inside the loop. Also, the 
apparent non-conservation of energy is explained, since the dispersive integral is 
done along the loop model imaginary part cut, where the initial energy is indeed 
sufficient to get on-shell constituents. In other words, the standard approach am- 
plitudes are constructed like in " old-fashioned perturbation theory" . In the context 
of quantum field theory, their natural framework is dispersion theory. 

Second, the spin projections are treated covariantly. The projector often seen 
in the literature is an approximation. In those approaches, the projector (1 -I- 70) P 
comes from the spin wavefunction of the positronium, while here it comes directly 
from the loop propagators. This is more adequate since it is the moving electron- 
positron entering the scattering process that must be constrained to the required 
spin state. The projectors appearing in (15) for particles in motion will introduce 
some new corrections to the positronium decay rate of the order of the binding 
energy 7^, i.e. (here "new" means relatively to the result found starting from 
(||) or any equivalent formula). 

Further, the covariant projectors have the important property of preserving 
gauge invariance. Indeed, when using (1 -I- 70) P as projector, gauge invariance is 
valid only if the constituents are taken at rest, i.e. k = (the so-called static limit). 
In other words, (15) is a gauge invariant positronium decay amplitude while (H) is 
not. 

Another consequence is that since the loop model appears at the root of stan- 
dard approaches, interesting alternative calculational tools can be built from it. In 
particular, it provides for an easy procedure to integrate the wavefunction exactly, 
keeping tracks of binding energy effects even at lowest order. This is the subject of 
section 3. 
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2.1.2 A comment about Bethe-Salpeter 

Viewing the decay process through a loop model is far from new. After all, it is the 
approach followed using Bethe-Salpeter (BS) analyses, where one starts with 



with r the off-shell scattering amplitude (given in (11) for n = 2). The loop 
structure emerges from the BS identification of the wavefunction as 



with Tbs the BS vertex function. All this is standard textbook material (see for 



As is well-known, the BS equation cannot be solved exactly, and one has to rely 
on some approximations. The most popular one for positronium is that of Barbieri- 
Remiddi , which can be viewed as a four-dimensional generalization of the usual 
Schrodinger wavefunction. To make contact with the three-dimensional integral 
representation for the decay amplitude (|^), the wavefunction is approximated as 
^' (?) - (5 (qo) ip (q) in (^ (see for example |ll|, (l5)). 

The novelty of our approach is to use dispersion relations to study the BS loop. 
This has great advantages over approximate methods, because our treatment is 
exact. Also, all the previous comments show the insights gained doing the reduction 
our way. In fact, starting with a four-dimensional BS form factor or vertex function 
Tbs ~ ^s(9o,q), the dispersion relations alone enforce qo = 0, i.e. the energy- 
dependence is set to zero. 

Another particularity of our approach is the treatment of spin, since we took 
the BS vertex Tbs with its spin structure replaced by 75 or ^ {FB{qo,fi) is a 
scalar), a replacement dictated by the bound state behavior under parity and charge 
conjugation. This allows us to identify the correct covariant spin projectors for the 
e'^e~ pair, and to preserve manifest gauge invariance. 

In conclusion, the present dispersion technique could help clarify some problems 
concerning the reduction of four-dimensional BS wavefunctions to three-dimensional 
ones (note that one should not confuse the present reduction of the BS wavefunction 
with the 4_D 3D reduction of the BS equation itself, extensively discussed in the 
literature, see for example and references quoted there). 

2.2 Analytical Properties of Positronium Decay Amplitudes 

We have just seen that the standard positronium decay amplitudes are given a 
natural interpretation in terms of dispersion relation integrals. Let us see the im- 
plications for orthopositronium decay into three photons. 

The loop model amplitude consists of six diagrams (figure 3) 




(16) 




example p^). 



A4 {o-Ps — > 777) 




(17) 



where the off-shell scattering amplitudes are 



r^''" (e+(iP - g)e-(iP + g) ^ 7 (^i) 7 (^2) 7 (^3)) 

-ze3(r^7 + rg,^^-frgr) 



(18) 
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with 
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+7^ 



and P = ^1 + ^2 + ^3- To compute these loop ampUtudes, we start by extracting 
the imaginary part, and we discover a completely different picture than the one 
described earlier for the two-photon parapositronium decay. 

Since the discussion is rather involved, the following summary of what we intend 
to show may be useful: 

1. There are six vertical cut diagrams V (figure 4a), with the properties: 

- reproduce the standard convolution-type decay amplitude; 

- are separately gauge invariant; 

- are badly behaved in the soft-photon limit, i.e. violation of analj^ticity. 

2. There are additional contributions arising from the oblique cuts V 
(figure 4b), with the properties: 

- are not of the factorized convolution-type; 

- are not gauge invariant by themselves, except for a constant Fb; 

- also badly behaved in the IR; the combination V + T> also fails to 

vanish in the SPL, except for a constant Fb- 

3. In addition to the loop model amplitudes V and T>, structure- dependent 
contributions S must be considered (figure 4c), with the properties: 

- are non factorizable, requiring a new form factor for the vertex Ps e+ 



We will spend some time to justify and explain these results, reviewing each 
point in turn. The strategy is to use dispersion relations as a classification scheme 
for contributions to the decay process, and fundamental principles (gauge invariancc 
and soft-photon limits) to characterize and constraint each type of contributions. 

Let us state our results more precisely, introducing some notations. We write 
the decay amplitude 



M {o-Ps (M^) ^ 777) = {M^^P (loop model) + M^"'^ (structure)) £^^£^£3, 



e 7; 



- act to restore gauge invariance; 

- interfere with both V and T> to ensure the correct soft-photon limits, 
(if Fb is constant, structure terms need not be considered) 

4. The various contributions should scale as 



^~0(a^), |~0(a^) 



= v+v+s 



Using the optical theorem, each contribution is expressed as 



V 




(19a) 
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= -tr 



i=i 

with the vertical cuts imaginary part 



ImV (P2) = el^el^el^j d^2M {o-Ps (P^) ^ e+e-) x M^^,';, (e+e 

(20) 

The integration is over the two-body e+e phase-space (i$2, and summation over 
spins is understood. ObUque and structure-dependent imaginary parts are 

(21) 

Im5i (P2) = el^el^el^j d^^M^j,.^, {o-Ps (P^) ^ 6+6-7;,) xA^^^^^^a (e+e" - 7^.7/3 

(22) 

and similarly for photon 2,3. The oblique cut amplitudes Mib are interpreted as 
bremsstrahlung amplitudes, hence the subscript. 

Gauge invariancc is expressed through Ward identities. Let us concentrate on 
the photon li for definiteness. First note that scattering amplitudes are all gauge 
invariant 

because the e+e^ are on-shell. This implies gauge invariance for the vertical cuts 

Zi^^V^-p = 

On the other hand, IB and SD amplitudes cannot be separately gauge invariant, 
since h,n.MiB-i 7^ 0- This serves as a constraint for SD amplitudes: we impose 

Provided this last equality holds, the positronium decay amplitude is gauge invari- 
ant. 

To analyze the soft-photon behaviors, let us expand each amplitude around 
/? = 0: 

n 

Selection rules alone sufHce to cancel IR divergences, i.e. vi = di = 0, and 

M {o-Ps (M^) ^ 7ii7i2 7(3) ' ~° Constant 
The implementation of the correct SPL requires in addition 

V2 + d2 + 82 = 
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which guarantees 

Similar behaviors are expected for photon 2,3. 

Let us now prove those assertions, and discuss their interpretations. 

2.2.1 Factorization properties 

By an immediate apphcation of the resuh of section 2.1, one can easily show that 
vertical cuts alone reproduce the standard three-dimensional convolution-type am- 
plitude: 

xTr y, (m- /(f') V^'^'P (e+ (fc') e" (fc) ^ 777) (to+ ^)] (23) 

with on-shcU e^e~. Note that it is important to treat spin projectors covariantly 
to preserve gauge invariance (if (m+ /(f) ^ {m— is approximated as (1 + 70) 
gauge invariance is lost, except in the static limit). 

Bremsstrahlung amplitudes generate a factorization problem for oblique cuts V. 
In fact, the emission of a photon by either the electron or the positron breaks the 
"symmetry" of the diagram (by symmetry is meant a configuration in which each 
lepton is carrying half the positronium momentum). Let us recall the implications 
of this symmetry for the vertical cuts. The first term in (|20| ) is 

M {o-Ps (P2) ^ e+ (f - q) e- (f + q)) ^ {q ~ iP, q + iP) (24) 

The energy conservation at that vertex and the on-shell condition for the e+e^ (in 
the o-Ps center-of-mass frame) amount to 

Siqo)s{\q\-^/pyI^) 

This allows us to write Fb very simply as Fb = C<j)QT {c^) (q^ +7'^) (12). Such 
simplifications do not occur for the oblique cuts, the energy-dependence of the form 
factor is not set to zero but to a complicated q-dependent value. We refer to the 
calculation of oblique cuts for p- Dm —f e^e~7 (section 4) for an explicit illustration 
of the problem. 

2.2.2 Gauge Invariance and Low's theorem 

Low's theorem states that an amplitude involving only neutral external bosons must 
vanish in the soft photon limit, i.e. when anyone photon's energy is taken to zero. 

This theorem must apply to positronium decay amplitudes. One must then 
conclude that standard amplitudes do not exhibit the correct analytical behavior, 
since individual scattering amplitudes contained in F (e^ (k) , (— k) — s- 777) ex- 
plode in the SPL. One can view this simply by noting the presence in F of electron 
propagators of the form 

i 

and when the ith photon energy Z° is vanishing, the propagator explodes because 
the constituents are on-shell. In fact, the situation is improved by a consequence of 
selection rules: when the six vertical cut amplitudes are considered, these IR diver- 
gences cancel and M (o-Ps 777) is finite. To understand this, remember that an 
IR divergence in (e+e^ 111) would be cancelled by IR divergences in radia- 
tive corrections to M (e+e~ 77), which cannot contribute to M {o-Ps 111)- 
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Nevertheless, and contrary to what is often beheved, selection rules do not suffice 
to enforce analyticity for the standard decay amplitude, because M (o-Ps — > 777) 
must vanish and not simply be finite. 

Looking back at the loop model, equivalent to the Bethe-Salpeter loop, a possible 
solution immediately emerges. Since only off-shell fermions can circulate in the 
loop, no analytical problem should occur. Taking into account all the cuts, both 
the oblique and the vertical ones, the correct SPL should be restored. In fact, this 
is wrong, because the loop model amplitude is not gauge invariant for o-Ps decays, 
and Low's theorem heavily relies on gauge invariance, i.e. Ward identities. Let us 
discuss this issue of gauge invariance more precisely. 

The source of the problem is in the momentum-dependence of the form factor. 
One can show that the loop model amplitude is gauge invariant if and only if 
the form factor accommodates linear shifts Fg (q) — Fg (q + h), and similarly for 
hih- For a constant form factor, the loop is gauge invariant, and no problem of 
SPL occurs (the loop amplitude is then obtained from the light-by-light scattering 
amplitude, for which it is well-known that gauge invariance guarantees correct SPL 

Equivalently, we can trace the problem down to the bremsstrahlung amplitudes. 
To visualize the situation, let us change the momentum parametrization and draw 
amplitudes for o-Ps e+ (pi) (^2) 7 ('1), as shown in figure 5, 

M'is;! = ieFB{p,+h,p2)l^u{pi)^^^^rviP2)^e^{P) 

MeFs {pi,P2 + h) l^u{p,) j-zMzljly (p) (25) 

Obviously, when contracted by li, the two amplitudes fail to cancel each other, due 
to the different momentum dependences of Fb : 

Iht^^iB;! = [Fb (pi + li,P2) - Fb {pi,P2 + h)] {u (pi) 7"« (P2)}e„ {P) + 

Thus gauge invariance is violated by the bremsstrahlung parts of the oblique cuts. 

Dispersion theory permits the precise identification of the bremsstrahlung pro- 
cesses as responsible for both analyticity and gauge invariance violations. The idea 
is now to add new structure-dependent (SD) contributions, such that gauge in- 
variance is restored. What we will show is that these new SD contributions also 
automatically restore the analytical properties for the whole amplitude. 

It could seem that the loop model has not achieved much. Worse, by considering 
it, we have not solved the analytical problems, and we have brought in new gauge 
invariance problems (since vertical cuts are separately gauge invariant). One could 
therefore think that it would be much easier to try to find new structure-dependent 



gauge invariant contributions to cure the standard amplitudes ( 23 ) , and forget about 
loops. 

Our point of view is that one must consider the loop model as the starting point, 
being strictly equivalent to the Bethe-Salpeter approaches (^6|). Another argument 
is that oblique cuts do a great deal of the work: for a constant form factor, they 
succeed at curing vertical cut analytical problems. The new SD contributions are 
therefore accounting for the non-constancy of the form factor. This is far more satis- 
factory. The loop represents the decay dynamics, and "nearly" succeeds at restoring 
SPL. The remaining violation comes from the bound state dynamics (non-constant 
form factor), and is cured by new bound state structure-dependent processes. 

Another advantage of the loop approach is that pseudoscalar para-state decays 
do not necessitate structure-dependent terms. As we will see, thanks to the ap- 
pearance of a pseudoscalar coupling 75, the loop amplitude is gauge invariant, and 
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well-behaved in the SPL. We will explicitly check that for the decay paradimuonium 
{fj,'^ fi~) into 76"*" e~, correct SPL is restored through the interference of oblique cuts 
with the standard approach vertical cuts amplitudes. This analysis is done in the 
last section, along with the analogous decay Ks — > je~^e~ . 



2.2.3 Structure-dependent terms 

The violation of gauge invariance by bremsstrahlung-type amplitude is a very well- 
known phenomenon in pion and kaon physics. One remembers that for a decay 
process like tt"*" e'^i^j, the so-called Inner-Bremsstrahlung (IB) amplitudes 

combine with Structure Dependent (SD) amplitudes to restore gauge invariance (see 
for example [E8|, or [pO| for a modern Chiral Perturbation Theory point of view, see 



also 1 34 1 for measurements of structure dependent form factors, and recent literature 
cited there). 

What we here propose is a similar mechanism for positronium. On general 
grounds, it was foreseeable that some interplay between the bound state dynamics 
and the decay process should appear, both ultimately involving photons. The present 
model allows us to give some constraints for the general form of those structure 
terms, exactly like for pseudoscalar pion decay. However, to obtain their exact 
forms is a very difhcult problem. 

To analyze the restoration of both the correct SPL behavior and gauge invariance 
thanks to SD terms, we invoke Low's theorem in its original form (see |l^, p8|): 

For a scalar decay process A'^ (pi) (^2) 7 (k) with both structure 

and bremsstrahlung contributions, the complete amplitude behaves in 
the SPL as 

M'^iA^^e+e-j) = f^e+e- (pi,P2) (Pi) ^^IpJ^^ (P2)| (26) 



if the Ward identity k^M^ = holds. Three important consequences: 

1- There is no constant term, i.e. O ^(fc*^)*^^, arising from the form 
factor. 

2- The form factor appears as evaluated at the same point. 

3- When the process involves only neutral bosons, A4 ^ O (fc") . 

To understand the content of the first point, let us expand Fb {pi +li,P2) in 
(25) around ll — 0. The expansion of the form factor will bring a new constant 



term 

FB{Pl+h,P2) —z ; Fb(Pi,P2)— ; 

2pi ■ li 2pi ■ li 

+2arf^^^|^+O0?) 
opi 

for some a. It is this additional constant term that must be cancelled by SD 
amplitudes to preserve gauge invariance, hence leading to Low's theorem (^6|). 

The second point can be rephrased as follows: near = 0, the complete ampli- 
tude behaves exactly as if the form factor was a constant (i.e. it is evaluated at the 
symmetric point (|2^ ) for all contributions). This is the key. The result is neatly 
expressed as 
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In the soft-photon limits, the loop amplitudes behave exactly like the 
constant form factor loop model, because structure- dependent amplitudes 
interfere with oblique cuts to cancel their dependences on the form fac- 
tor variations. As a result, positronium decay amplitudes vanish in the 
SPL, since for a constant form factor loop model, oblique cuts interfere 
destructively with vertical cuts. 

This result is quite powerful, because the dynamics of each contribution is man- 
ifest. This opens the way to an evaluation of their respective scalings in a and 7, 
as we now discuss. 



2.2.4 Scaling and NRQED 

At first sight, all contributions (vertical, oblique and structure) should be of com- 
parable numerical size, all of them being of the same order in a. This is in fact not 
at all the case. A very special feature of electromagnetic (loosely) bound states is 
that the binding energy can be expressed in terms of a 

7 w ma/2 

It is one of the great virtue of NRQED to provide a consistent framework to re-order 
the perturbation series according to powers of a"7'" ck"+™. 

Another virtue of NRQED is that the starting point is the Bethe-Salpeter loop 
model (^6|), i.e. the same loop model as (|l7|). This would mean that all the cuts 
should somehow be contained in NRQED calculations. If it is true that lowest 
order NRQED starts by setting go = in (|l^) (because kinetic energy is subleading 
compared to momentum), systematic corrections to this approximation are then 
calculated, and the four-dimensional loop is "reconstructed". In other words, the 
lowest order NRQED result is our vertical cut results, while oblique cuts correspond 
to higher order corrections. 

Our approach could then be viewed as complementary to NRQED perturbative 
expansions. Dispersion relations put forward basic principles, and should provide 
interesting constraints on the systcmatics of NRQED. 

Structure-dependent contributions are not contained in NRQED. Their form is 
unknown and so is their numerical importance. Their existence means that the issue 
of gauge invariance for the NRQED expansion is irrelevant, because hopeless. This 
means that the " reconstruction" of the four-dimensional BS loop is not constrained 
by gauge invariance. However, by showing that structure-terms are contributing 
at the order a^, one proves that the spurious gauge dependence and violation of 
analyticity of NRQED is also of that order, hence validating the present O (a^) 
theoretical predictions for the rates. 

Let us review the different contributions. What follows is not a rigorous demon- 
stration, but rather a qualitative analysis. Nevertheless, the exact evaluation should 
not alter much our conclusions. 

i- The vertical cuts alone are sufficient to reproduce the known lowest order 
decay rate 



r (p-Ps — + 777) — 



2 (7r2 - 9) 



This result is obtained in the static limit, i.e. with Fb oc 5'^'^^ (q). Let us also give 



the differential rate Olll 



dV {o-Ps 777) 2a^m 



dxi Qtt 



2(2-a:i) 2(l-3:i)^i 



In(l-xi) 
(27) 
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with xi = l\/2M the reduced photon energy. The spectrum ( p7| ) iUustrates the 
violation of analyticity. Near zero, it vanishes linearly 

dT {o-Ps 777) x^^o 5 
d^, 3^^^ 

while we would expect a xf behavior instead. Indeed, the phase space alone is 

dr (o-Ps -f-fj) 2a^m 



dx 97r 



[2xi] (29) 



Since the amplitude should behave as xi near xi = 0, the differential rate should 
exhibit a x\ behavior in the SPL. As we have seen, the additional contributions 
restore the correct SPL. 

ii- The oblique cuts are suppressed by three effects: they disappear in the 
limit of vanishing binding energy (7 — > 0) , the dispersion relation integration range 
is reduced compared to the one for vertical cuts, and the three-photon phase-space 
enhances vertical cuts against oblique ones. 

The first point is obvious: if the e+e^ pair emerges on-shell from the bound 
state, there is no room for oblique cuts. In fact, this limit is pathological, because 
oblique cuts are still there, but contribute only when one of the photon energies is 
strictly zero, to enforce SPL. Anyway, their complete contribution should somehow 
be proportional to the binding energy 7, i.e. the mass default of the e^e~ inside 
the positronium. 

To explain the second point, consider the dispersion integral (19b). The lower 
bound Smin(/i) is increasing as increases. Indeed, if the photon is carrying away 
some energy, the initial bound state mass squared s must be greater to create on- 
shell intermediate e+e~. Therefore, oblique cuts should be decreasing functions of 
the photon energies. 

The first two effects combine to suppress oblique cut contributions everywhere 
except near zero photon energy. Therefore, oblique cut contributions to the decay 
rate will be further (substantially) suppressed by the phase-space (|^, peaked for 
high photon energies. 

Due to all those effects, one could state that oblique cuts should contribute at 
the order 7^ « a^, i.e. T) « a^V . Since the argument is qualitative, it could happen 
that they are enhanced for some reason (or further suppressed), but we do not see 
any compelling evidence of such pathologies. 

iii- Finally, structure terms can be estimated by taking the conservative view 
that their importance is just what is required to preserve gauge invariance and 
SPL. Since they compete with oblique cuts, they are of the same magnitude, i.e. 
O (a^). However, structure terms account for the variation of the form factor. For 
a Schrodinger-like form factor 



(q2+72)2 



The variations will be a much more peaked function of q , leading to a smaller 



result for the dispersion integral (19c). All in all, we should gain a power of the 
binding energy S « 7I?. In terms of the fine-structure constant, structure term 
should be of order , i.e. S ~ a^V. 

Of course, this is not very convincing. If the pion decay can be used as a guide, 
structure terms should indeed be highly suppressed. However, the very small mass 
gap between the parapositronium and orthopositronium could lead to huge enhance- 
ments, if one imagines structure term dynamics as driven by intermediate virtual 
parapositronium. 
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In conclusion, we argue that, provided that no pathological effects come into 
play, oblique cuts could be considered as small O (a^) corrections, and structure 
dependent terms neglected at the present theoretical precision level. Anyway, no 
answer to the orthopositronium lifetime puzzle could be viewed as entirely satisfac- 
tory as long as no definite estimation or calculation of structure-dependent processes 
has been achieved. 



3 Form Factor and Parapositronium Decay 

We have established the correspondence between (10) and (15). Let us construct 
an alternative, but equivalent, dispersion procedure specific to the two-photon case 
that will be used in explicit calculations. 

By computing the trace in (10), the tensor structure factorizes and we are left 
with 

M {p-Ps 77) = Sme'e^'''P"lij2,ael^eia (M^) (30) 
and the decay rate is expressed simply as 

r {p-Ps ^ 77) = mna^m'^M^ \X (M^) |^ (31) 

In this equation, X (M^) can be viewed as an effective form factor, modeled as the 
electron-positron loop with the coupling Fb- There is only one term in X (M^) 
since the direct and crossed amplitudes are equal under q —^ —q, i.e. an allowed 
variable change as Fb {q^,P ■ q) — Fb {q^, —P ■ q) , and we write 

1 1 1 



(27r)" (q ~\P) - m? {q + ^Py - {q ~ \P + h 



(32) 

It is to evaluate the effective form factor X (P^) that we will now use dispersion 
techniques. Remark that the factorization of the tensor part is interesting, since 
gauge invariance is manifest, and that X {P^) is convergent while the amplitude 
(10) is superficially divergent. 

The factor 77 is introduced because there is a subtlety in the above factorization. 
Indeed, there is an arbitrariness in the choice of variable for the dispersion integral. 
This situation is well-known for the photon vacuum polarization: 

(fc2) ^ (fc2^M'> _ ni (fc2) = (^g^- - U2 (fc2) (33) 

where one writes a dispersion relation for Hi (fc^), which is less divergent than 
112 (fc^) due to the factorization of the tensor structure. The analogue of (33) for 
our case is 

M {p-Ps ^ 77) = Sme'e^^-'P^hJ^^^el^etli {M') = Sme^e^^'^'^^^et^eLIa {m' 

We will choose to write a dispersion relation for X2 (M^), and this corresponds to 
the choice rj = P'^/M'^. This seemingly arbitrary choice is in fact necessary to 
be consistent with the general development, see (p^, (p^. In other words, the 
dispersion relation was built on the whole amplitude, so it is clear that the photon 
momenta appearing in the tensor structure were also incorporated, i.e. they were 
reduced as h^p h^p x \/P^/2. That is the reason why we must include a factor 
P^ into the effective form factor X (P^) ■ 
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Let us give a general expression for X (-P^) as a dispersion integral. The loop 
integral X {P^) has an imaginary part obtained by cutting the propagators 



1 p2 

lniX(p2) = iil- 



2TriS 



Fb 



-\P) 



{q-\P + h) -m^ 



Proceeding exactly like in the general case (see appendix), using an unsubstracted 
DR, we reach the two equivalent forms 



X (M^) = 



C<t>o 1 
2Af 2 327^2 

Ccf)o 



+ 00 



dsJ^(s/4-?7i^)ln 



1 



1 - 



1- x/1- 



4m^ 



2AP 



d^q ^, 2\ 



(34a) 

(34b) 

q| cos 9q 

The first expression is 



where we have set Fb = C(j)oJ- (q^) (q'^ +7^) (see |i; 
obtained by integration over the angular variable 9q. 

This is our third representation for the same decay amplitude : the first is the 
loop integral (10), the second is the well-known amplitude (15) (or (||), with no 
approximation for projectors) viewed as a dispersion integral for the amplitude, 
and the third is the present dispersion integral (34a,b) for the effective loop form 
factor X (M^). All three procedures are strictly equivalent to each other. 



3.1 Decay Rate in the Static Limit 

Let us first analyze the static limit, i.e. 7^ ^ for the form factor. To compute the 
decay rate in that limit, we do not need to specify Fb yet. We just need to know 
that the function J- (k^) is normalized to unity and behaves as a delta function of 
the momentum in the limit of vanishing binding energy: 



(k2) = 1, lim T (k2) = (2^)3 5(3) (j^) 



(35) 



By setting !F (q^) = (27r)'^ S^^'' (q) in (32), we must recover exactly the lowest order 
decay rate 

r {p-Ps — > 77) — -a^m (36) 

In some sense, the static limit serves as a boundary condition for the behavior of 
J". 

Using (34b), we get 



X (M^) 



C(f)o 

2AP 



1 



|q| cos 6' 



C(j)o 

2AP 



Importantly, this result is independent of the binding energy : the loop do not 
introduce any corrections in the static limit. The decay rate in that limit is therefore: 



1 ^ m'^ 
T{p-Ps^jj)^-a^mi—C' 



(37) 



with |0o| = a^m^/Sn. It remains to match C such that purely kinematic correc- 
tions vanish {M factors in the above formula arise from products like (li ■ I2) — 
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M"^ /2 and from the 1/2M decay width factor, while m comes from electron propa- 
gators in the loop and from the wavefunction 0o). With the definition C = \fM jm, 
the decay rate is exactly ^a^m as it should. In other words, the value for C ob- 
tained by matching (15) and (||) is such that no correction arises from factors M/m 
in the static limit. 

To conclude, let us repeat that we have not specified the form factor. This 
means that any form factor which has a three-dimensional delta function limit for 
7^ ^ gives the correct lowest order decay rate ^a^m. In the following, we shall 
present three forms, all built on the Schrodinger momentum wavefunction. 



3.2 Schrodinger Form Factors and Binding Energy Correc- 
tions 



The formula (34a) is a very simple expression: the decay rate (31) for any form factor 
is obtained by a one-dimensional integral. This is a very nice calculational tool 
that we will now use. We will go through different calculations of F {p-Ps 77), 
obtained for specific choices of Fb (or equivalently. (q^))- Namely: 

^n{e) = (38b) 

(q^ + 7 ) 

Tiii(q^) = ^arctan— x Ti (q^) (38c) 

|q| 7 

where 7^ = m? — Af^/4 w related to the binding energy Eb — M ~ 

2m = —ma^ /A. All these form factors satisfy the delta limit property (35), but are 



differently peaked near |q| = (see figure 7). As a result, we expect that all of them 
will reproduce the lowest order decay rate ( |36[ ) , but will introduce some 7-dependent 
corrections proportional to their spreading around |q| = 0. Such corrections can 
then be expressed in terms of the fine-structure constant, yielding a corrected lowest 
order rate of 

r [p-Ps ^ 77) = ia^m (1 + Sr (T)) 
The physics of each form factor is discussed below. For now, let us use the for- 



mula (34a) to get the effective form factors and decay rates corrections 5r {^i.ii.iii) 



For the first one, we find 



Ij (M^) = —p^- arctan — (39) 



The integral needed for this calculation is (see [21]) 
dx 



In 



1 + Vi^ 



and its derivative d/dxo- The decay rate is 



"'"^^^ 2 arctan^ , (40) 
y/xo - 1 



1 f Am?\ [2 
F/ {p-Ps 77) = -a-'m —r - arctan — (41) 



2 V M"^ ) V"' 27 



where we used |0o| = ct'nt' j^-K and C = \J M jm. By expanding this result around 
7 = 0, and expressing corrections as a series in a, we recover the standard result as 
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zeroth order 

Ti{p-Ps^-f-f) = la^m ( 1 - - + la^ + O (a^)^ 

2 \ TT 8 / 

« ia^m(l- 0.637a + 0.351a2 + O(a^)) 

Proceeding similarly for the other two, using ( ^0| ) for JF//, and numerical integration 
for J- III, wc find 

Tii{p-Ps~^-fj) « ia^m (1 -0.25a2) 

Till {p-Ps -> 77) w ia^m (l - 1.73a^) 

If we naively combine the present corrections with radiative corrections up to order 
In a: 

r,.p. = ^^l-(^5--j- + 2a^ln- + .p(^)j 
the theoretical value is modified as 



Lowest-order 


Corrections 


Decay Rate, with 


Form factor 




Radiative Corrections and 5r 


Static Limit 





7.9895 X lO'-'sec-^ 


^/ 


-0.637a 


7.9527 X lO^sec-i 


^// 


-0.25a2 


7.9894 X lO^sec-i 


^iii 


-1.73a2 


7.9887 X lO^sec-i 


Experiments 




(7.9909 ±0.0017) x lOysec"^ 



where the experimental measurement (nl) is also included for comparison. 



3.2.1 Discussion 

In view of the above table, it appears that the first form factor JF/ leads to a too 
small value for T {p-Ps 77), since it introduces new order a corrections. The 
problem is in the form factor, which does not converge fast enough towards the 
static limit delta. In other words, for a given 7^, J^i (q^) is not enough peaked 
around |q| = 0. On the other hand, the more peaked J^n and J^m form factors 
lead to acceptable corrections. 

It is time now to consider the physical content of each form factor. In summary, 
the J^i is just the Shrodinger momentum wavefunction for the bound state, the J-'n 
form factor has no clear signification and only serves the purpose of illustration, 
while J^iii can be viewed as a Coulomb binding corrected Schrodinger form factor, 
in the spirit of the Sommerfeld factor. Of course, one could argue that the form 
factor to be used must be linked somehow to the Bethe-Salpeter wavefunction. In 
fact, it is the first form factor Ti that emerges from the reduction of the Barbieri- 
Remiddi wavefunction through dispersion relations (15). Despite of that, we will 
argue that it is the J-m form factor that must be considered if one is interested in 
O (a) accuracy (this discussion is much inspired from p| ). 

The fact that the exact integration of the !Fi form factor leads to excessive 
order a corrections i5r {^i) is not new. Let us stress first that those corrections 
are 7-dependent (see (41)), i.e. it is a binding energy effect. Such effects, being 
non-perturbative, are quite difScult to be dealt with. Even in the O (a) calculation 
of ||ll[, much space is devoted to the discussion of such effects. Double-countings 
plague bound state calculations, and it is quite possible that the present excessive 
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O (a) corrections should be discarded because they are already taken into account 
in O (a) calculations . 

Indeed, the standard approach to these excessive corrections is simply to discard 
them. Either they are assumed to be part of higher order radiative corrections to 
the scattering process (jl^, or alternatively, they are considered as part of the 

relativistic corrections to the bound state wavefunction |Q, iQ. In any case, they 
are interpreted as accounting for some modification of the intermediate e'*"e~ state: 
either as emitted from the bound state, either as entering the scattering process. 

This dual view is further enhanced by the following: those additional corrections 
disappear in the static limit, i.e. when 7^ —^ 0, so they appe ar a s originating from 
the bound state dynamics. Alternatively, looking back at ( |34lj ) , it appears that 
the electron propagator between the two photons plays a crucial role in introducing 
the additional corrections. If it is frozen (i.e. replaced by 1/m), the corrections 
disappear. In other words, one could view the decay process instead of the bound 
state dynamics as responsible for the excessive corrections. 

From the point of view of the BS equation, things are quite clear. The Schrodinger 
momentum wavefunction is built as the Coulomb photon exchange ladder approx- 
imation (see figure 8a). Relativistic corrections to the wavefunction must be con- 
sidered to account for the covariant photon exchange. This corresponds to the 
relativistic corrections point of view. The radiative corrections view arises by con- 
sidering the binding graph (figure 8b). As shown for example in this O (a) 
radiative correction for the scattering amplitude reproduces the lowest order decay 
amplitude, because its Coulomb photon exchange part has already been accounted 
for in the construction of the wavefunction. 

What we propose is to view the third form factor as a representation of the 
modification of the intermediate e^e^ state. When dispersion relations are used, 
one consider intermediate states as asymptotic states, and integrate over the corre- 
sponding phase-space. This is the application of the optical theorem for absorptive 
parts. Let us assume that, due to the modification of the asymptotic Hilbert space 
by the long-range Coulomb interactions, a correction factor must be introduced in 

(p_P5 ^ 77) ^ J d^q (q2) xSiq)xM {e^e+^ ^ 77) (42) 



M 

and that this factor is 



S{q)^p- arctan M (43) 

|q| 7 

In some sense, one could think of S (q) as the Sommerfeld factor of Harris and 
Brown [8]. Due to long range Coulomb interactions, the wavefunction at contact 
(f)o is " renormalized" by 

by using v » a, i.e. |q| >> 7, with v the center-of-mass velocity. Correspondingly, 
S (q) for small binding energy is 

^(^)-|qr2. 

The factor of 2 arises because S (q) is defined at the amplitude level, while the 
Sommerfeld factor is at the decay rate level. 

The same arctangent factor arises in the work of 0, fill on O (a) radiative 
corrections, by a detailed analysis of the binding graph and its Coulombic part 
(which is, after all, an initial state Coulombic interaction). 
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In conclusion, if one is interested in the lowest approximation (i.e. the ma^ 12) 
the static Hmit is just fine. If one wants the lowest approximation to O (a) accuracy, 
one must consider it as arising from the O (a) binding graph, as demonstrated by 
Adkins p] . In other words, by taking only O {a) diagrams for the scattering 
ampHtude e+e~ — s- 77, one finds 



1 c / ^ a I _ TT 



2^ 



r {v-Ps 77) = ^rna' ^1 - - 1^5 - — j j + O [a ^ 

The 1 in the above expression is the lowest approximation to O (a) accuracy, arising 
from the O (a) amplitude. To reproduce it using our effective form factor method, 
the Coulomb corrected Schrodinger form factor J-ju is to be used, as in ||ll|. 

What we gained using our method is a better understanding of its link with 
the Sommerfeld factor, i.e. asymptotic Coulomb interactions. J-m being more 
peaked than J-j, the suppression of the intermediate phase-space is interpreted 
as a manifestation of the iterative structure of Bethe-Salpeter construction of the 
wavefunction. 

Therefore, the method of is nothing but the rephrasing of that of Harris 
and Brown [||. However, the method of ||ll| or the present integration of is 
surely more appropriate, because one does not rely on static limits, hence avoids 
the well-known static divergence 1/v of radiative corrections to e~e^ 77 p6[ |. 

3.2.2 Improving the perturbation series 

It is interesting to reverse the previous argument. Instead of curing the form factor 
so that one avoids double-counting, it would be more adequate to modify radiative 
corrections, keeping the J-'i form factor. Indeed, the form factor accounts for most 
of the O (a) corrections. We found 



Yj [p-Ps [1-2^ + {a^)) « ^ (1 - 0.637a + O {a^)) 

while complete corrections are 



p-Ps 



1 _ ( 5 - ^) - + O (a^) ) « (1 _ 0.806a + O (a^)) 



2 V V 4/77 
In other words, the decay rate can be expressed as 



a^m 4m^ / 2 M 



arctan— 1- 3 ] - + O (a 



a 



2 NP \t: 'i'l J \ V 4 / TT 

a^mAm^ (2 mV , / 9nn 

■arctan— (l - 0.170a + O (a^)) (44) 



p-Ps 



2 M2 \tt 27 J 

This last form is very interesting because binding energy corrections (i.e. 7-dependent) 
are singled out, while the resulting purely radiative corrections are much reduced. 
This means that one could, at least in principle, express the decay rate as non- 
perturbative binding energy corrections times a rapidly converging perturbation 
series of radiative corrections. This is exactly what we have achieved to order a. 

It is worthwhile to note that ( |4^ ) also accounts for most of the O (a^) correc- 
tions. Omitting logarithmic corrections, the decay rate is given by 

a 7Ti 

Tp_p, w (1 - 0.806a -H0.52a2 + O(a^)) 

am 4m /2^^^^^^M\ „ q 170a -I- O.OGa^ O (a^)) (45) 



2 Af2 \t: 27 
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This appears as coincidental (especially in view of the many interferences occurring 
among O (a^) corrections, see but it is nevertheless intriguing. It may be 

the signal that an improved basis can be built for the positronium perturbative 
calculation, if one avoids NRQED systematic a expansion of binding energy effects. 
(Remark that this formulation does not cause any problem regarding logarithmic 
corrections; the modification amount only to a reduction of the coefficient Cp of the 
In a term to -3.919(1) from -7.919(1)) 

Of course, it would be very interesting to build the corresponding expression for 
orthopositronium. Unfortunately, as extensively discussed, both oblique cuts and 
(unknown) structure-dependent contributions render the calculation highly chal- 
lenging. 



4 Paradimuonium Decay into 76+6 

The present section concerns paradimuonium, the singlet /i^/i^ bound state ||27|| . 
This state has not been observed yet. The reason why we consider the decay p- 
Dm — > 7e+e^ is that it is the simplest process where a photon is kinematically 
allowed to have a vanishing energy (Note that the positronium decay p-Ps — > 71^77 
via a Z° has the same dynamics, but has not been observed either). Therefore, we 
will be able to illustrate that the loop model has a correct soft-photon behavior in 
the case of parapositronium or paradimuonium decay. 

The decay p-Dm — > e+e^7 is shown in figure 9. The simplicity of the pro- 
cess comes from the pseudoscalar nature of the paradimuonium, which allows a 
manifestly gauge invariant treatment throughout (hence the absence of structure- 
dependent terms). The amplitude for p-Dm — > e"'"e~7 can be directly obtained 
from the positronium to two photon decay amplitude (|3^), (32): 



M (p-Dm e+e-7) = 8me3e^''''^fcpi,< (fc) (P) (p')} ^ ^^2^ n ^^g) 

-\- le 

with X ~ 2kP jM the reduced photon energy, M the dimuonium mass and m the 
muon mass. The effective form factor loop integral is given by 

J (P x) — ri / —Fb 

^ 7 (27r)^ (g-iF)'-m2 (g+ip)'-m2 (g-ip + A:)'-m2 

where rj — P^ and — P^{1 ~ x) (the only difference with p-Ps 77 being 
that t'^ 7^ 0). Gauge invariance is ensured by the factorized tensor structure, i.e. by 
the antisymmetric Levi-Civita tensor. 

From the decay amplitude, the differential width is obtained as 

— ^ ^ = — — m^M^ |i(m2,x)| pix,a) (47) 

with the phase space spectrum 



p(x,a) = Jl-- [a + 2(l-a;)]- ^ (48) 

V i — a; [1 ~ X) 

where a = Am^/M^, and the bounds on x are [0, 1 — a]. This spectrum is highly 
peaked just below 1 — a, as shown on figure 10. 

We are going to calculate the integral I (M^, x) using dispersion relations. Be- 
cause P (i.e. X ^ 1), two different types of cuts will contribute (figure 11). 
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We have shown that considering the vertical cuts is strictly equivalent to the de- 
cay amplitude calculation done using formula (15) where F is now the (amputated) 



scattering amplitude /i^ (fc') /i (k) 77* 



"7 with on-shell muons. 



What we are going to show is that the amplitude for p-Dm — s- e+e"7 has a 
correct soft-photon limit only if we take all the cuts into account. By this we mean 
that each cut gives a contribution to the amplitude that behaves as a constant when 
the photon energy goes to zero. The combination of the vertical and oblique cuts, 
on the contrary, forces the amplitude to vanish in that limit. We thus recover the 
analytical behaviour expected from Low's theorem for the decay p-Dm — > 77* 
involving only neutral bosons. 



4.1 Dispersion Relations and Soft photon limit 

In this section, we will analyze properties of the effective form factor I (P^,a;). To 
keep the discussion as general as possible, we will not specify the form factor Fb- 
i- The absorptive part is found by cutting the relevant propagators as 



2 (27r)^ 

dS 
2(27r)^ 



■Fe 



2TriS((q-^P) -nA2TTiS((q+^P 



2m 



S (^(q + ip)^ - 2t:iS ((g ^ + kf - m^) 
(g-iP)'-m2 



for the vertical and oblique cuts, respectively. By a straightforward integration, the 
first expression gives (s = P^): 



Imli (s, x) = --^Fb (<7o = 0, |q| = ^s/A-mA In 



1 + yjl - Arn^/s 
1 - a/1 - 4m2/s 



while the second one cannot be completely integrated without specifying Fb'- 



{s - 4m^) 
(49) 



Iml2 (s, x) 



77 1 



S \&-KX 



dqo 

qo 



Fb qo, |q| = yqi + qoVs + s/4 - 



with the bounds given by 



ymin — 4 ' ^ 



1 - 



4m 



s{l-x) j ' * 



4 



1 - 



s{l-x) 



4m"' 
" 1-x 
(50) 



(51) 



Interestingly, the second cuts contribute for qo =/= 0. This is in sharp contrast with 
the two-photon decay, since there only the first cuts exist. Further, approximated 
bound state wavefunctions where a S (qo) appears cannot be used (see [11] to p^), 
and one should revert to the full four-dimensional Bethe-Salpeter wavefunction (for 
example the Barbieri-Remiddi one |24[). 

ii- In the soft photon Hmit, the combination ImXi (s, a::)+ImX2 (s, x) behaves 
as a constant when x despite the fact that each cut diverges in that limit (see 
the 1/a; in ( ^9|) and (50) ). This will guarantee that the imaginary part of the whole 
amplitude vanishes in the soft photon limit thanks to the presence of kp in the 
tensor structure of (^). The following simple formula gives the behaviour of the 
imaginary part when a; ^ 



ImXi (s, x) -I- ImX2 (s, x) -ttt- 

S iDTT 



fE 



\/s — Am"^ 



1 - 



0{x) 
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where 



/b = FB(0,q2 = s/4-m2) 



f's 



QFb (go,|q[ 



+ 



9i=B(?o,|q[ 



+0,|q|->-^s/4-m2 



To get a feeling of the underlying physics, note that for a constant form factor 
Fb = the absorptive parts are 



ImXi (s, x) 



S IQlTX 



In 



1 



1 



Iml2 (s, a;) = — 



s 167ra; 



1-^1 
1+^/1- 



4m^ 



e{s- Aw?) 



In 



s{l-x) 



1 



s(l-x) _ 



Arrr 



The relative —1 factor is responsible for the destructive interference when x 
Explicitly, the SPL is 



0. 



Imli (s, x) + Inil2 (s, x) 



m Z?Const 



s 167r ^1 - 4mVs 



ill- The dispersion integrals are 



x] = - 



ds 



1 /-"^ ds 



M2 



ImXi (s, a;) 
ImX2 (s, x) 



(52) 
(53) 



We expect that the dispersive part will inherit the SPL properties of the absorptive 
parts. 

The oblique cuts are seen to be decreasing functions of x, and disappear when 
^ Q [x ^ 1). Further, in that limit Ii {s,x) ^{s), with X(s) the two- 
photon effective form factor (32). At the decay rate level, the oblique cuts are 
really negligible because the phase-space spectrum p {x, a) is highly peaked around 
X = 1 — a. In fact, independently of the form of Ti {s,x), provided it is not too 
badly behaved, we can set X (Af^, x) I\ (M^, l) in (47) to get an estimation for 
the ratio as 



(IT [p-Diii — c+c--) /dx 
r {p-Dm 77) 



a 



37r |J(M2,1)|- 



-p{x,a) « —p{x,a) 



(54) 



with a = Am\/M'^. This is a rather universal result: the highly peaked spectrum 

wipes out the details of the dynamics, and the same ratio holds for other pseu- 
doscalar (7r° e+e~7. Kg e+e~7,...). Integrated over x, this gives 



nth 



r {p-Dm 



-7) 



r {p-Dm 77) 



(55) 



while for example Rj^o ~ f- (15.3) ~ 0.012, in good agreement with experiments 
[34]. 

If wc took some time to recall known stuff, it is to emphasize that the present 
discussion of p-Dm e~^e~j is really to be taken as an illustration of some suppres- 
sion effects. For the more interesting orthopositronium decay, similar kinematical 
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effects are expected. However, comparing the respective phase-space photon spec- 
trums (figs. 6, 10), the suppression of obhque cuts is seen to be less effective in 
orthopositronium. 

We will now evaluate the contribution of each type of cuts for the case of the 
Schrodinger momentum wavefunction form factor (38a). In doing so, we will gain 
a better understanding of the analyticity restoration at the level of the differential 
decay rate for p-Dm — > e+e~7. 



4.2 The vertical cuts reproduce standard approach results 

To precisely show what happens when the oblique cuts are forgotten, we compute 
here the rate with the vertical cuts only. This is the result one would obtain starting 
with (15). 

From the first cut imaginary part (^9|), the decay amplitude is found through 
an unsubstracted dispersion relation. The calculation is very similar to that of the 
p-Ps — > 77 decay amplitude, and we get 



M3 X 



2 M' 
— arctan — 
TT 27 



(56) 



Inserting the result for Ti into the expression of the differential rate with C' 
M/m' and |0of 



a^m^ /8tt, we get: 



dT (p-Dm 



7e+e ) 



dx 



a^m / Am? \ 
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M 


6tt V M2 j 




arctan — 


TT 
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Note that, independently of the value for 7, the photon spectrum is always linear 
in X for a; w 0, due to the 1/x in (56). This is an incorrect soft photon behaviour 
since Low's theorem requires a cubic spectrum instead (the amplitude behaves as x, 
and an additional x comes from phase-space). The 7-dependent corrections cancel 
in the ratio 

dT (p-Dm —* je^e^) /dx a p{x,a) 



r {p-Dm — !■ 77) Stt x'^ 

(compare with (54)). For completeness, the total rate is simply: 



(57) 



r [p-Dm 



7e+e ) 



r {p-Dm 
Numerically, 



77) 



a 
3^ 



-Vl - a (a - 4) + 21n 



1 + VT^ 
1- VT~^ 



r {p-Dm 



76+6 



r {p-Dm — > 77) 



(18.7) 



to be compared with (55). 

4.3 The oblique cut contribution to the rate 

We have just seen that the vertical cuts suffice to reproduce the lowest order decay 
rate. It is interesting to investigate how the oblique cuts can restore the analytical 
behaviour of the spectrum without affecting this lowest order evaluation. 

Unfortunately, as mentioned earlier, the form factor to be used to compute 
oblique cuts imaginary parts (and thereby their contributions to the decay rate) is 
not the simple Schrodinger momentum wavefunction, because the energy-dependent 
part is not set to zero. To keep the discussion as concise as possible, we will 
nevertheless use this Schrodinger wavefunction, i.e. 



Fb (90, q) = C(t)o 



q2 ^ ^2 
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We could say that the standard approach neglect of subleading kinetic-energy de- 
pendences is taken. In fact, since our purpose is only illustrative, the present ap- 
proximation is not very important. The cancellation leading to the correct SPL can 
still be observed. Of course, the exact value of the oblique cut contributions to the 
decay rate will not be exact, but this is not a urgent question for now (the decay to 
e"'"e~7 is about 1, 5% of the rate to 77 and the oblique cut corrections are at most 
of a few percent, while p- Dm has not even been observed yet!). 

To be clear, basic principles (gauge invariance and analyticity) require the intro- 
duction of oblique cuts, but practically, for low order evaluations, their exact form 
is not relevant. On the contrary, a precise definition of the form factor is necessary 
to tackle seriously the calculation of o-Ps 777, but this is left for future work 
(because it may be more useful to uncover the remnants of oblique cuts in standard 
NRQED expansions first). 

When using the Schrodinger form factor, the imaginary part of the effective form 
factor has the constant SPL 

ImX, (.,.) -MmX. (.,.) Ti^^p:!-^ (33) 

We will now integrate the dispersion relation (53). Obviously, this integral is quite 
complicated and not very interesting for the present purpose. Instead, we revert to 
numerical evaluation of 2'2 (M^,a;) as a function of x, and compare it to (|56|). 

In the figure 12, we plot the vertical cuts contribution (dashed line), normalized 
to get a constant as 

— / ,.9 x fxAP \ ^ ,,,, X 2 M 
I, (m2, x) = — — Ii (A^^ x) = - arctan — 

The complete effective form factor, with the same normalization (solid line) is 

T{M^x) ^ i^^yi,(M^x)+J,{M^x)] 

2 M 1 xM^ /■+°° ds , ^ , , 
■ arctan - — 1 / —r Iml2 (s, x) 



TT 27 TT CSo Jiinl S- M'^ 

l~-x 

The plots are drawn for m = 1,M = 1.9975 (i.e. 7 « 0.05) and for m = 1, M = 
1.99975 (7 w 0.016). Note that the physical value 7 « ma/2 corresponds to 7 ~ 
0.0036. For x ^ 0, one can verify that T (M^, x) 0. 

The figures show that near zero the two cuts interfere destructively in order to 
maintain a correct analytical behaviour for the whole amplitude. Away from x — 0, 
the oblique cut contributions are strongly suppressed relatively to the vertical one, 
and this suppression increases as 7 decreases. As can be seen on the graph, it is 
typically for x < 7/771 that the oblique cut contributes. 

Therefore, we can summarize by giving a simple representation of the different 
contributions to the amplitude. From the figure 12, one can see that the behaviors of 
the vertical cuts, the oblique cuts and their combination are quite precisely modelled 
as 

Ii - -, X2 , ' ,,,, ^ X = Xi ' - ' 



x{x + -f/M) \x + j/M 

As a consequence, the spectrum behaves as 

dT{p-Dm^ -fe+e~) 2 , . z( 1 V ..^s 
- |X| p(a;,a)-x'' — ■ — (59) 



dx I I A-v w \x + -flM 
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i.e. a linear spectrum when 7 — > 0, and a spectrum when x is smaU. The effect 
of 7 7^ is therefore to soften the photon spectrum and sUghtly reduce the total 
width. This behaviour is exactly what we postulated in a previous work p9[ |. 

In conclusion, the oblique cuts have a small contribution to the decay rate com- 
paratively to the vertical cuts. However, their presence is essential to guarantee the 
analytical properties of the amplitude expected from Low's theorem. 

4.4 Overview of Kg — > e+e~7 

The restoration of analyticity by a cancellation among the vertical and oblique cuts 
is presented in the context of the kaon decay —> e"^e~7 via a charged pion loop 



(see 28 to and also [^). In this analysis, the form factor describing the 
K —^ TTTT vertex is taken as a constant. This decay is interesting since there is a 
close similarity between this hadronic decay process and the present QED bound 
state decay. The main conclusion is that even if charged particles are present 
in intermediate states, the amplitude has to vanish in the soft photon limit, as 
expected from the fact that the decays Kg —>■ 77* or p-Dm 77* involve only 
neutral external bosons. Further, this Kg decay provides a physically sensible 
process where one can analyze the constant form factor assumption, since for such 
a loosely bound system as the p-Dm a constant form factor cannot be realistic. 

Note that the modern approach to this decay channel is Chiral Perturbation 
Theory, see for example ||2^ , |Q . 

4.4.1 The pion loop model 

Details can be found in ||l9|, so let us simply review the main steps. The decay 
amplitude is modeled as a charged pion loop, to which photons are attached either by 
the one-photon or the two-photon seagull coupling (see figure 13). Gauge invariance 
is thereby enforced, but is manifest only after the loop integration has been done 
using dimensional regularization. At this stage, the amplitude reaches a form similar 
to (PI) 

M{Kl^je+e-) ^ (x" ^ ^ ^ (M^ x) x 

(47r) 

(fc) {" [P) IvV [P )} 

with t=p + p' = P — k, rriTr the pion mass, x the reduced photon energy 2k^/M 
and M the kaon mass. J-'{s,x) is the Feynman parameter integral 

J^is,x)^ f dy f ' dz- ^. ^"'^^^ , (60) 

Jo ^ Jo 1^4{a-b)zy + 4by{y~l) + ie 

with the definitions a = s/Am^^h = s{l — x)/Am'^ (in the K rest-frame). This 
Feynman parameter integral plays exactly the same role as X (M^, a;) , i.e. that of an 
effective form factor (the only difference being that the momentum loop integration 
has been done, and the imaginary part is obtained using the prescription ie, instead 
of by cutting propagators, this is anyway strictly equivalent). The differential rate 

dV [Kl ^ e+e-7) Idx 1 ^|2 , , 

r (Ag n+TT ) Stt'^ a^V 1 - a^r 

with Qtt = 4mJ/M^, Oe — Aml/AP, me the electron mass. The phase-space is given 
in (48) with the same bounds x G [0, 1 — Ce]. 
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Using the same pion loop model, the two-photon decay rate is similarly obtained. 



and 



r {Kl ^ 77) 



1 



\T(M\l)\' 



r {K° TT+TT') 7r2 o2yr^ 

As previously mentioned (Q), the ratio of the two electromagnetic modes is ap- 
proximately given by 

dr(K°s^e+e--f) /dx a\T(M^,x)\'' , , a , , 
r(i^o^77) 3^ |^(Af2,i)|2 3tt^^ 

Anyway, as we will now discuss, the precise x behavior is important regarding the 
enforcement of the correct SPL. 



4.4.2 Dispersion Relations 

The effective form factor T (Af 2 ^ x) is calculated using dispersion relations |^l[] , |^ . 
Its imaginary part consists of the vertical cuts (see figure 13) 



IraTi (s, x) 



27rm2 2m^ 
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4m^ 



{l-x)dl 



is - Am"^ 



and the oblique cuts 



Im (s, x) = 



27rm2 2to2 



■In 



' s(l-z) 



a-x) ji 



' s(l-x) ^ 

(compare with ImX(s,a;)). Using an unsubstracted dispersion relation 
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J^{M^,x) = Rejr(M2,; 
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Im (s, a;) 



(61) 



for A'f2 < 4m2 . The final expression for T can be obtained by analytical continua- 
tion for any value of Af 2 



T{M^,x) = - 
with 

fix) 



1 



1 



1 



5(2;) = 



2(a-6) ' ^^_,fy2^fi-)-f(b)) + H9ia)-9{b))] (62) 



arcsin2 (^/x) < a; < 1 

— (In + \/x — 1) — ii-Tr) x>\ 

arcsin (^/x) < x < 1 

(In (^ + Va:- 1) - \iT^) x>l 



This result agrees with that of O (p**) Chiral Perturbation Theory [29], [|3^, or with 
a direct evaluation of the Fcynman parameter integral (|60|). 

4.4.3 Soft-Photon limits 

Individual contributions of each cut to the imaginary part are divergent as ~ l/a;2. 
However, their combination is finite in the soft photon limit a; ^ 0. Indeed, applying 
L'Hospital's rule twice, we get 



Im jFi (s, x) + ImJF2 (s, x) ■^=" 



1 



s2 ^1 - 4m2/s 



(63) 
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This constant SPL is not altered by the dispersion integrals, and the effective form 
factor behaves as a constant for very low x 

(note that for < 1, the imaginary part of (64) is the same as in (|6^)). Since there 
is a factor k in the tensor structure, the whole amplitude vanishes in the SPL, as 
expected from Low's theorem. 

Let us turn to the soft photon behaviour of the differential width. Since the 
amplitude behave as x for low x, the resulting spectrum is in (x^ from the squared 
amplitude and a x from phase space). In other words, if we had forgotten one cut, 
the resulting spectrum behaviour would have been divergent as 1 /x near zero instead 
of vanishing like x^ . Compared to the paradimuonium and orthopositronium cases, 
it is even more crucial not to forget oblique cuts, because of the IR divergence that 
would be spuriously generated. 



4.4.4 Non-constant Form factor 

For a non-constant amplitude M. {K^ tt+tt^) = F (^{q + fc)^ , [q — t)^^, structure- 

dependent terms must be supplemented to enforce gauge invariance [g8|. By an 
analysis strictly parallel to the orthopositronium one, the origin of the problem can 
be traced down to the bremsstrahlung amplitudes 

{K% ^ ^+ (pi) vr- (p2) 7 {k)) = F ({p, + kf ,pl) ^-F (pj, {p^ + kf) ^ 

The effect of the structure terms is to restore both gauge invariance and analyticity 
by accounting for the variations of F. 

In conclusion, the pseudoscalar decay K'g — > e^e^7 is seen as an interesting 
system, being the low energy mesonic analogue of dimuonimn decay to e^e~7. 
The dispersive techniques introduced to deal with electromagnetic bound states 
are the standard tools used to described e~^e~'y. Also, the same general 

electrodynamical principles must hold. Namely, gauge invariance and analyticity 
in the SPL are shown to be implemented by the same mechanism, i.e. interference 
between imaginary part contributions. For non-constant form factors, structure 
terms come into play to account for the variation of the form factor, to preserve the 
cancellation among contributing cuts in the SPL. 



5 Conclusions 

In this paper, we have used dispersion techniques to analyze positronium decay 
amplitudes. In this way, emphasis is put on basic principles, i.e. gauge invariance 
and soft photon limits. This provides interesting new insights into bound state decay 
dynamics. Indeed, basic principles implementation requires the non-perturbative 
treatment of some binding energy effects. What our work has shown is that even 
if NRQED scaling arguments are a very powerful calculational tool, the inherent 
breaking of covariance hinders the manifestation of basic principles. In other words, 
binding energy effects are usually perturbatively expanded along with radiative 
corrections in the NRQED approach. The present work can therefore be understood 
as complementary to NRQED, and may emphasize its underlying model-dependent 
assumptions. 
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Some consequences are worth repeating here. First, no approximations were 
needed to reach standard three-dimensional convolution-type decay amplitudes (ex- 
cept of course the neglect of oblique cuts and structure-dependent amplitudes). 
This provides a more concise (and consistent) scheme to deal with spin structures 
and energy-dependences of the Bethe-Salpeter wavefunction. As a by-product, we 
showed that some formulas usually quoted in the literature are approximations, 
missing some of those O (a^) corrections they are meant to evaluate. 

Also, our method emphasizes the fact that the exact integration of lowest order 
decay amplitudes contains O (a) corrections. This is not an artifact of our model. 
Rather, this is a theoretical issue; a great deal of work in positronium decay calcu- 
lation being devoted to avoid double-counting. In our view, the most appropriate 
formulation is to consider the corrections arising at the lowest order as binding en- 
ergy effects, and therefore to express them in terms of 7 (see (H])). Further, this 
reordering of the non-relativistic perturbation series may be much more adequate 
to deal with QCD bound state like quarkonia, for which little is known about the 
wavefunction. 

Concerning orthopositronium, a consistent picture of the decay process is built. 
Gauge invariance and analyticity are at last correctly implemented. Further, by 
exploiting the fact that some destructive interferences must occur among the various 
contributions, bounds can be extracted on their magnitudes. Unfortunately, only 
lower bounds are given, but one could reasonably expect that those bounds serve as 
estimates, and conclude that the present O (a^) NRQED calculation is complete. 
However, one should keep in mind that no proof of this conjecture exists at present. 

Let us make a comment about this optimistic conclusion. We have shown that 
the decay process at lowest order in a requires the introduction of both subleading 
contributions (oblique cuts) and additional structure-dependent contributions. This 
requirement is to be understood at the level of basic principles. In particular, 
gauge invariance is valid only if structure terms are taken into account, i.e. gauge 
invariance does link photons inside the bound state to radiating photons (emitted in 
the decay). We have seen that NRQED proceeds by a perturbative reconstruction 
of the Bethe-Salpeter loop. If gauge invariance is taken as a constraint in this 
reconstruction, some contributions may be missed, because the loop is ultimately 
not gauge invariant. Of course, since structure terms are invoked to restore gauge 
invariance, and since they should be of order O (7'^) ~ O [a^) , such questions 
are not important practically at the present level of precision. Anyway, it would 
be interesting to perform a thorough analysis of gauge-dependent contributions 
arising in NRQED, in order to get an independent estimate of structure-dependent 
contributions. 

The central achievement of this work is to single out and characterize the 
structure-dependent contributions as responsible for a breakdown of standard NRQED 
at the O (a^) level, or even sooner in pathological cases. By breakdown is meant 
that considering only a two-body e'^e'^ Bcthe-Salpeter wavefunction is no longer 
sufficient (structure-dependent terms could in principle be obtained from a two-body 
to three-body Green function) . One should understand positronium calculations as 
a two step process: the wavefunction is obtained from BS analyses, while the decay 
process is calculated using NRQED. So it is really the Bethe-Salpeter basis that is 
at stake, not the non-relativistic effective theory. In our view, theoretical advances 
now require going beyond the two-body Bethe-Salpeter approach. No answer to 
the orthopositronium lifetime puzzle could be given before the completion of such 
progresses. 

Acknowledgments: C. S. and S. T. acknowledge financial supports from 
FNRS (Belgium). 
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6 Appendix 



This appendix contains the demonstration of the assertion that 



Dispersion relations constructed on the vertical cut imaginary parts of 
the loop model reproduce standard convolution-type amplitudes. 

We will particularize the discussion to the parapositronium decay into 27, for 
which there are vertical cut contributions only. Let us emphasize that the whole dis- 
cussion of this section is readily extended to any para- or orthopositronium vertical 
cut contributions to the decay amplitude. 

We first compute the imaginary part of (10) for an arbitrary initial squared mass 
P^. Considering the two possible cuts (figure 2), we obtain ImT by replacing the 
two propagators on each side of F^'' by delta functions 



ImT (P^ 



ImX^'^(p-Ps^27)£t^£L 




After a straightforward integration over and |q|, with P — we reach 
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In the course of the derivation, the delta functions forced q'^ ~ and |q| — 

(65) 



^ P"^ /A — m? . In other words, the electron momenta are 

2 02 




— ,±q with UP±q] = — -h\' 



This kinematics is to be understood in the trace evaluation. The angular dependence 
arises from the relative orientations of q and the photon momentum li. Note also 
that the relation (65) cannot be satisfied for the physical value P^ = < Avn? . 
This is obvious since the loop cannot have an imaginary part for the physical bound 
states, its constituents being always off-shell. From the kinematics ( |65| ) one can 
prove that the factors on both sides of T^^'^ are true projectors, which serve to 
enforce gauge invariance in the expression 



t 
2 



Indeed, those two projectors play exactly the same role as external spinors when 
demonstrating Ward identities. 

The real part will now be calculated using an unsubstracted dispersion relation 



RcT(Af^) = 
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+00 



ds 



A/2 



ImT is 
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(66) 



where it is understood that P^ should be replaced by s everywhere, i.e. scalar 
products that will appear when evaluating the trace should be expressed with the 
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kinematics defined for an initial energy s. Since Af^ < Am?, the principal part can 
be omitted and T (M^) = ReT(M^). Now let us write the form factor in the 
general form 

Fb = C<i>oT (q2) (q2 + 72) = c^oT (s/4 - m^) • (s - M^) /4 (67) 

with 7^ = — M"^ /A and 0o the bound state wavefunction at zero separation. 
Then (66) can be written as 



r(M2) = C<j>o I ^ ds [ '^Tis/A-m 

Jim^ J AtT 



647r2 

2^"'J'^ \^)« ^ + J ^1^1^21^ 

Let us transform the s integral back into a |q| integral, keeping in mind the 
constraints obtained when extracting the imaginary part. Using = s/4 — m^, 
ds = 8 |q| c?|q|, the decay amplitude dispersion integral is 

where, as the notation suggests, it is understood that any appearing in the am- 
plitude must be replaced by 4|q|^ + 4m2. In particular, ^/P^~(q} can be replaced 



by 2i?q with Eq — \/|q| +m'^. This amounts to consider the scattering ampli- 



tude with incoming on-shell electron-positron having momenta (^P (q) ± g)^ = 
(since q'^ — Q). Note the fact that Eq > Af/2, apparently the energy is not con- 
served. This is not surprising since the present formula is a dispersion integral, done 
along the cut where (q) > 4to2 Finally, in view of the kinematics, we introduce 
k = ^P (q) + q and k' = ^P (q) - q (hence Ek = Ek' = Eq and k = -k' = q) to 
write the amplitude simply as 

^ (^^') = T / TT^sTTT ['t^oT (k^)] Tr {75 (- ^ + ^) T'^'^ (fc, k'M)i]^ + m)} el^el, 
^ J (27r) iii/k 

(68) 

where T^'^ {k,k' is the amplitude for on-shell e (fc) e"*" (fc') scattering into 27. 
Gauge invariance is present due to the two projectors, well defined since k^ — k'^ — 
777.2. rpjj-g gj^jjg Q^j, demonstration, and T (M2) = M. {p-Ps — > 77). 
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Figure 1 : The loop model direct and crossed diagram for the decay p-Ps ^ '-fy 




Figure 2 : The two cuts contributing to the imaginary part 




Figure 3 : The loop model diagrams F^s'' for the decay oPs — > 777 
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Figure 5 : Bremsstrahlung processes and structure-dependent contribution 
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Figure 6 : Photon energy spectrum for the decay o-Ps — > 777. 




Figure 7 : Form factors q^^(q^), as functions of with 7 = 0.5 
{Ti plain, Tu dashed and Tui small dashed lines) 
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Figure 8 (a): BS vertex equation (Coulomb photon ladder approximation) 
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Figure 8 (b): The binding graph. 
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Figure 9 : The loop model amplitudes for p-Dm — > e+e 7 
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Figure 10 : The photon energy spectrum p{x, a) for a = 0.01 
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Figure 11 : Vertical and oblique cuts imaginary parts 
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Figure 12 : The vertical cuts (dashed line) and vertical plus oblique cuts 
(plain line) contributions, for 7 0.05 and 7 0.016, as functions of the 
photon energy x. 
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